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ABSTRACT 


At the present time the primary method of obtaining solutions to 
non-linear differential equations is by means of the digital computer 
and numerical techniques. A method is here proposed to find an approxi- 
mate mathematical expression through the use of Laplace Transform tech- 
niques. Thus, the Laplace Transform concept is extended to the solution 


of non-linear differential equations. 
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CHAPTER 1 


1.0 Introduction 

Laplace Transforms have long been employed in the solution of 
linear differential equations. Baycura [ref. 1] proposed a method of 
obtaining an approximate solution to a non-linear differential equa- 
tion by using a linearized transform to obtain the approximation. In 
a later paper [ref. 2] it was proposed that an exact solution might be 
found by developing an expression for the Laplace Transform of a non- 
linear term. The result was that the solution was not exact, but very 
Similar in form. 

These results are examined in detail] here, and a general expres- 
sion for the non-linear transform is developed. Further, the concept 


is extended to non-homogeneous equations. 


1.1 Applications 


A method of obtaining a good approximation to the many non-linear 
differential equations which occur in all engineering fields would be 
of great utility to the engineer. Having a mathematical expression 
which closely approximates the solution would enable the engineer to 
analyze the system, and to optimize the design more rapidly, and 


without having to rely totally on digital computers. 


CHAPTER 2 


2.0 Development of the Laplace Transform for Non-Linear Functions 


In order to use Laplace Transform techniques in obtaining an 
approximate solution to a non-linear differential equation, expres- 


sions for the transform must be derived. 


2.1 The Type of Equation Defined 


Consider a general non-linear differential equation of the form: 
(eye eta Mines ment se) (2.1.1) 


where x) (4) is the men derivative of x(t); k indicates the derivative 


raised to the kt" power; and x(t) is x(t) raised to the ne power. 


2.2 The Method of Successive Integration by Parts 


To derive an expression for the transform, use is made of succes- 
sive integration by parts. Consider two time dependent functions, 
x(t) and y(t). The integral of their product is given by: 

fxy dt=x fy dt-x ff y dt dt +x fff y dt dt dt - ... 


(-1)PxhM) i y dt|" dt ( 2a) 


where ] y i is an n-fold iterated integral. This equation can 


now be used to find a series expression for the Laplace Transform of 


x(t) by letting y(t) = = This: 
[xe"s* dt & x oe ah | les" dt dt + x | [fest ase ake wins Soe 
0 O O O 


Wie) 


After integration the series becomes: 

co -st, <-6t7  :- Seite 
[xe"s* ee . Xe | _ xe xe 

S S 

O O 








2 =f. (2.2.3) 








e) 


In evaluating the integral at the upper limit the condition must be 


such that: 


Limit ye St = 0 (DD. A.) 


t > @ 


and similarly for all derivatives of x(t). Thus, evaluating the 


integral at the lower limit, the result is: 











[xerst fe oe x(0) x(0) tives beee 39) 


0 
and since: 


e [xes* dt (2B) 
O 


by definition, an infinite series expression for the transform has 


been derived. 


2.3 Derivation of an Expression for the Transform of Derivatives 





To derive an expression for the transform of a derivative of 
x(t), the method just illustrated is used again, with x(t) in place 
of x(t). Thus: 


oO 


[est dt = {erst “lino x{ fers 


OO 


eS* dt dt 4 i ife** de dt dig— --- 


oe 











O © 
XO) ek) FeO ee (aah, 
( oe ee 


Now, by adding and subtracting x(0), and then multiplying by s/s, the 











result is: 
[xs dt = 50) + SU] + A + | - x(0) (2.502) 
O 


The term in brackets is just X(s), so: 


A[x(t)] = sX(s) - x(0) (ers) 
a result which is well known. In a Similar manner it is found that: 
Bae) |S SAS eee) 0) (2.3.4) 
This method can be continued and it soon becomes apparent that the 
th 


expression for the transform of the m” derivative is: 


ol Ex™ (t)] = s™(s) - sM™x(0) - sM-2K(0) - 06 - gx'M2) 9) 


_ yi™1) (9) (2.3.5) 


2.4 Transform Expression for Integral Powers of a Function 


The functions considered thus far have all been linear; now the 
functions which cause the non-linearities will be considered. First, 


consider the function squared. 


[xzerst ‘alia = [x(xe"*)¢t 


O @) 
2 x |xe"S* dt - if fxs diedie. ain st dt dt dt - -:: 
° ° (2.4.1) 
Now, let 
gis xx" dt (2.4.2) 
0 


The integral is the same as Equation (2.2.6), so: 


I, mex OFS) (Qe 3 } 
Now let 
FR i| [rest dt dt (265) 


and if the inner integral is examined: 


[xe"s* dt = xfens* dt - i{ [erst dt dt + x{{{-$* dt dt dt - -:- 
-St *-St -—St 


= - - - 4-4 - (2.4.5) 





This result is then put in Equation (2.4.4) to obtain: 


s2 
O O 0 


Ip = (Est dt + ts die 4 byes" a = (2.4.6) 


Each of the integrals can be evaluated by using the results of Section 





2235. S05 
* X(s ] ] 
In = (of : bs 1,]54s) - (0) + S]eeKls - sx(0) - io) 
(2.4.7) 
If this is expanded and like terms are collected, the result is: 
_ x(0)X(s) 
aaa (2.4.8) 
Continuing in this manner it is found that: 

LE? (t)] = [x(0) + MDs Os TGs) (2.4.9) 
But the term in brackets is just sX(s), so: 

A x et) ee sXe (S9) (2.4.10) 
Thus, the transform of the square of a function has been derived. Now 
consider: 

[xte"s* dit= x|xze"st ae = if arent dt dt+--- (2.4.11) 
O O O 


The first integral has already been found, so, from Equation (2.4.10): 


x] xe st dt = x(Oysxots) (2.48.72) 
O 
The second integral is evaluated in the same manner as Equations 
(2.4.4) through (2.4.8). Thus: 
. i [x2e"S* dt dt = x40) SKS) (72.8 13) 
O 


The remaining integrals are evaluated in a like manner and the outcome 


1S: 








at be (tidus sx?(s)[*{0) , (0), x(0) , | (2.4.14) 


Once again, the term in brackets is recognized as being sX(s), so the 
expression can be written as: 

of [x?(t)] = s?X*{s) (2.4.15) 
If higher powers of the function are examined it soon becomes apparent 
that: 

Kix" (t)] = sl xs) (2.4.16) 


2.5 Transform Expression for Powers of a Derivative 
Since many equations involve a derivative which is raised to a 
power, it will be useful to find an expression for the transform of 
such a function. First consider the square of the first derivative: 
[(x)2e"S* dt = ike" st dt - i{ fre" Stat i (2.5.1) 


O O O 
The integral in the first term has been evaluated in Section 2.3, so 


this term becomes: 


i fiers ake eG GUS Se ate) (255R2 ) 
O 


The second term is evaluated in the same manner as Equations (2.4.4) 


and (2.4.13). Therefore: 


CO 


: X| fest dt dt = U%rsx(s) - x(0)] (2.5.3) 
O 
The remaining terms in the series are evaluated in the same way, so the 


transform becomes: 


{L(t F} = [sx(s) - «(0)2 fx(0) + 0) , X10) , J 








wary (sia (0 (s) [24 OM | (2.5.4) 


S S2 








But the second bracketed term is equivalent to sX(s) - x(0), so the 


transform finally becomes: 


F{lx(t)}4} = sEsx(s) - x(0)]? (2.5.5) 


10 


It should be noticed that this last expression may be written as: 

L{La(e) I} = sfCb(t)}? (2.5.6) 
If the foregoing method is applied to the third power of the first 
derivative, the result is: 

Lftx(e3} = se{feecey} (2.5.7) 
If higher powers of the first derivative are examined, it can be 
Shown that: 

xfricey a} = s Upper} (2.5.8) 
Now powers of the second derivative will be examined. First consider 
the second derivative squared: 

[(a%e"s* dt = ifrerst dt - x{fiers* dt dt + ++. (2.5.9) 

O fe) fe) 
The integral in the first term is known to be s*X(s) - sx(0) - x(0). 


The second integral is treated in the same manner as Equations (2.4.4), 
(2.4.13), and (2.5.3) and becomes (1/s)[s*X(s) - sx(®) = x(0)]]). “iis 


it is found that: 


. -tegy (4) | 
ff} = Gi(o) + MO 4 XD) y ..yLs2x(5) - sx(0) - X(0)3 


S S 
(2.550) 





The first term in brackets can be written as: 


HO. 0. | > eo , (0), X(0) | Tr 





S 2 Ss Se 
(0) 4 UO) gs «(0)| ~ sx(0) 
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wn 
nN 
< 
— 
=) 
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ae 





S 
= 5°40) + x(0) | = §°x(0) = sien) 


s{s? 2) + x10) A | - sx(0) - x(0)} (2.5.11) 


This is recognized as being s[{s*X(s) - sx(0) - x(0)]. Thus: 


L{ir)13 = sfet Lx (t)]} ? (2e5mi2) 








a 


2.6 Summary of the Transforms 
If Equations (2.3.5), (2.4.16), (2.5.5), (2.5.7), and (2.5.12) are 
examined, it can be seen that a general expression can be written for 


the transform. Thus: 


{E06 (eI = seh ey} (2.6.1) 


wheremim=n0, | ,92eeSreeeees andens= lue2, 35-4 
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CHAPTER 3 


3.0 Application of the Transforms 


The use of the Laplace Transform expression just derived will now 


be demonstrated by obtaining solutions to several equations. 


3.1 Application to a Non-Linear, Homogeneous Equation 


Consider a homogeneous, non-linear differential equation of the 
form: 

e(t) + Ae(t) + Be?(t) = 0; e(0) = V (3.1.1) 
Equation (2.6.1) is now used, and the transformed equation is: 

sE(s) - e(0) + AE(s) + BsE*(s) = 0 (3.1.2) 
This expression is now rearranged to give: 

BsE7(s) + (s + AJE(s) - V=0 (3 ee) 
This equation is quadratic in E(s), so it may be solved for E(s) 


_ through the use of the familiar quadratic formula. Thus: 


(sy) wc es ee (3.1.4) 


S 
Now this result is rewritten to put it into a more convenient and 


useful format. Thus: 


E(s) = el : + f i (in) ‘ues 


Now, if the positive radical is used and is expanded in a binomial 
series, the resulting expression for E(s) is: 


E( ) ey BV2s 4 Zev osS* 5BAy%s3 + (3 ] 6) 


Thus, an infinite series expression has been obtained for E(s). Now 


the inverse transform is to be found for each term of the series [ref. 3], 


iis: 


and the resulting expression is: 


~ eI 242 
e(t) = Ve ah pve : e| + (Bvt)? - ry, At | 


9/9 Onto AS tee fain 
= (BVt) 3 - “oe 8 - THe + - (Ble) 


An infinite series approximation has been found for e(t). However, for 








the original equation, an exact solution is known [ref. 4]: 
SIGE verMth - Me" = ih 43.128) 


By expanding the term in brackets, the two solutions may be compared. 


Thus, after using a binomial expansion, the exact solution is: 


2\72 
es vert h js we ; i : ole ™ ‘ 1}? 


3y73 
+ BA fooAt a)? a (3.1.9) 
Now, if each of the expotentials is expanded in a series, the result is: 


= \e7 At At wuhet2@ gbatdengue 
e(t) = Ve javelin = SE - + | 





2 6 
votveet(t = at + WE MEE) | 

(2h. 110) 
Notice the similarity between this equation and the approximate solu- 
tion in Equation (3.1.7). Each sub-series in the approximate solution 
is truncated, since the second term of the s-domain expression yields 
two time-domain terms; the third term of the s-domain expression gives 
three time-domain terms; and so on. 

From an analysis of the parameters A, B, and V, it is apparent 

that BV < A, and further, the smaller BV is with respect to A the 
better the approximation will be. Table 1 shows a comparison of the 


exact and approximate solutions for two sets of B, V, and A. Figure | 


14 


Shows a graph of the exact and approximate solutions when A = 1.0, 


B= "mlesemand ¥ = le. 


TABLE 1 
A. SAS TRO ees, VS 
Time beast Approximate Error 
0.0 1.000 1.000 0.000 
0.5 O07 0.473 0.034 
10 Ome2 0.8257 0.022 
eS 0.161 0.174 -0.013 
2,10 0.094 0.118 -0.024 
Cw 0.056 0.057 -0.001 
5. U 0.034 -0.010 0.044 
B.. A =-1,0, B= G.125, V Sale 
0.0 1.000 1.000 0.000 
Oi 0.578 0.577 0.001 
Ae 0.341 0.344 -0.003 
a5 0.203 0.214 -0.011 
2.0 Omle2 @ 1139 -0.017 
20 0.074 0.095 -0.021 
340 0.044 0.067 -0.023 
250 Oa027 0.048 -0.021 
4.0 0.016 Q.035 -0.019 
4.95 0.010 0.026 -0.016 
9. 0 0.006 0.019 -0.013 


{| dundL4 








Spuodas ‘] 
0°S “ey ars 0°2 O°L 
——— on | >) woo 
azewt xouddy iL 0 


O°L =A = (0)? 
Gcl 0 = dq 
OL= WY OF (3)29a + (2)av + (3)a 203 
#91 


| 
| 

| SUOLZNLOS azewLxXOUudde pue 4Dexd 3YyY 4O UOSLUedWO) 
EYE 
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3.2 Application to a Non-Linear, Non-Homogeneous Equation 


Consider an equation of the form: 

Mv + kv? = Mg (3u2.A,) 
This is the equation of motion of a falling object with the retarding 
effect of air resistance considered. It is assumed that the retarding 


force iS proportional to the velocity squared. If the equation is 


written as: 
v+ hv? -g=0 (69272) 
where 
v(0) = 0 (oe2is) 
and 
= k/M, (soc 


then an approximate solution to the equation can be found by applying 
Equation (2.6.1) as before. Thus: 

sV(s) - v(O) + hsV?(s) - 2 = 0 (372.5) 
This expression is rearranged to get: 

hs*V2(s) + s?V(s) - g = 0 (Sma2 sb) 
Once again the result is a quadratic equation in V(s), and it is 


solved for V(s) by using the quadratic formula. Thus: 


W(s) = gpia|- 8? # RY + aahs? | (3.2.7) 


This equation is now rearranged to get the same form as that used in 


Equation (3.1.5) 


a at + san) (302ce) 


If the term in parentheses is expanded in a binomial series and the 


positive sign in front of the radical is used, the result, after 


Simplification, 1s: 


“hos 20h | Scie seal ciclmee oie 
Vs) = Sp = + BED Say Mare Aa . G.2.9 


and the time-domain solution is obtained by finding the inverse trans- 


form of each term of the series. Thus: 


“ 243 3 2p 5 4) 347 S,449 
v(t) = gt - a + eon ik gt 14 ' I fe 20) 


So an approximate solution has been found. The exact solution is 


known [ref. 5]: 


qgi%il1 - o-2(gh) *t 
v(t) = 24-2 (3.2.lae 
2 
1 +e e(gh)"t 
where (g/h)% is the terminal velocity of the fallina object. If the 
L 
numerator and denominator of this expression are multiplied by elgh)?t 
then: 
oe 3)" tanh (gh)’*t (3.2.12) 


The hyperbolic tangent can be written as an infinite series, and in 


this form the exact solution is: 
Ss 


v(t) = (3) (any _ Chl C0 ae lig) -](3.2.13) 


To compare the two solutions, Equation (3.2.10) is rewritten as: 
vi 
By, 


ven # (Bon aytonyFer , (a) etonper (5 partgnie ,.. 
(3.2.14) 
The series are similar in form, but each term of the approximate solu- 
tion is smaller than the corresponding term in the exact solution. In 
order to get a numerical comparison, consider the falling object to be 
a spherical rock with density (op) of 2600 kilograms/meter’. The 
following parameters must be defined and evaluated [ref. 6]: 


M = pV 


Kk = 025 0. Go A 


d 


where M is the mass of the object, V is its volume, and A is its 


surface area. The other constants are: 


i 


1.293 kilograms/meter*® (Actually, pa is a function of air 
pressure and temperature, but it 
will be assumed to be very nearly 
constant for this analysis.) 


Pa 


Cy = 0.5, the draq coefficient for a sphere 


A = nD*, the surface area of a sphere 
3 
V = a the volume of a sphere 
g = 9.8 meters/second*, the acceleration due to gravity 


Thus, if the sphere has a diameter of 0.1 meter, the exact solution is: 
V(t) = seen4 tana 0.27% (3. 2s) 
and the approximate solution becomes 


(One 7s" 
6 


v(t) 2 36.4|0.27t ; (0.2782. (0.;2%6)’ » (Oper | 


"eo  Ooe ” 2aeeG - 
(3.2.06) 
The comparison of these two expressions is shown in Table II and a qraph 
in Figure 2. In the qraph, notice that the exact solution asymptoti- 
cally approaches the terminal velocity 36.4 meters/second, while the 
approximate solution overshoots that value, and after about seven 
seconds decreases monotonically. This is due to the fact that the 
highest power term in the series was preceded by a negative siqn, and 
as t became large this term overpowered the others. The first term 
omitted after the series was truncated was positive, so if it were 
included, the graph would end by increasing monotonically. The hump 


is caused by the lower degree terms being predominant when t is small. 


z. 
= 


GO PM 


Ge 


So ca S&S 


> > WwW 
Oo oOo © 
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TABLE II 


Exact Approximate Error 
0.000 0.000 0.000 
4.884 4.885 -0.001 
9.590 9.682 -0.086 
Tse e021 -0.182 
17.945 18.679 -0.734 
21.413 7 -1.318 
24.373 26.403 -2.030 
26.845 29.649 -2.804 
28.872 32.440 -3.568 
Seo 10 34.757 -4 247 
31.816 36.598 -4.782 
32.847 alipeekers' -5.12] 
33.656 38.879 -5.223 
34.287 So) BN -9 050 
34.776 39.328 -4,552 
35 3 38.788 -3.635 
35.444 37.574 -2.130 
35. 668 35.399 Oreos 
35.840 343750 4.090 
39.972 Zoe OG 10.204 
36.073 Wo. 077 e996 
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3.3 Application to an Equation with a Third Degree Non-Linearity 


Consider now a non-linear differential equation of the form: 

x(t) + wex(t) + hx?(t) = 0, x(0) = A, x(0) = 0 ceca 
After using the transform expression of Equation (2.6.1), the trans- 
formed equation is: 

So Xe) + (s*4w°)X(s) - As =0 (Gi. 3e2) 
This is in the so-called "normal" form of a cubic equation, in that 
there is no term involving X?7(s) present. Thus there are solutions 


X (s), X (s); and X (s); where: 


_ iB 
x = A+B, and x x = - &(A+B) + —5 (A-B) (.3.3) 
and where: 
| 
_\IA A? So 
A,B = \ons * \anzs% * “27has& (3.3.4) 


It can be seen that expanding this expression and determining which 
root to use is a formidable task. However, it is known from Equation 


(22,5) that: 


X(s) = A. + ait) + mo) pes hoes 5) 


x< 
om 

© 
~~” 








and if the initial conditions from Equation (3.3.1) are used and if it 
is assumed that all higher derivatives can be neglected at t = 0, then 
Equation (3.3.2) may be written as: 

hAsX?(s) + (s*+u))X(s) EeAs = 0 (3.3.6) 
Thus, the degree of the equation has been reduced by one and the solu- 


tion can be obtained in the same manner as in the previous two sections: 


a (Ss tuoi | 4hA?s*? }% 

X(s) Ee 1 +/1 + (s?Fue)? (393.7) 
Then, expanding the term in brackets as in previous solutions, the 
result is: 

ea hA’s° ZhoAgse 





sttue ~ (s@tuz)* * (sttuz)> 7 (3.3.8) 
22 


Now the inverse transform is found for each term of the series. For 


the first term, then: 
& | cad = A cos wt (302. 9) 
O 


The inverse transforms for the succeeding terms are more difficult to 
evaluate, as no expression of the form s'/(s?4w2)" can be found in the 
tables of Laplace Transforms available. However, one handbook [ref. 7] 


does have the transform: 


4 n-% 
(stxa8)r quent “nglvol! (3.3.10) 
O 2° "T(nju, 


Further, the table [ref. 8] has the relation that: 


1 


s"g(s) = f\) (4) + + p(k) (gyno Ik N Bly 25 038ecren( al. 


k=o 
er Fk) (Q) = 0 for k= 0, 1, ..., n-l. Thus, these two relations can 
be combined to yield the desired inverse transforms. So, for the 


second term of the series: 


S 
g(s) = 722, 2)3 (3.3.12) 
(s*+u)) 
and its inverse transform is given by: 
rat? J#( uot) 
2 P(3)u, 
Then, the Bessel Function is given by [ref. 9]: 
1 
7 9 \2(sin wot 
J3(ut) = ec a - COS Wet (3.3.14) 


Upon substituting this expression into Equation (3.3.13), the resultant 
inverse transform of the second term is, after applying Equation 


CORSO 


— 
Y 
w 


Rs : 242 
f(t) = (3H t sinw t tut cos wt) Goose ls) 


oe 
of 
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This method is then applied to the other terms of the series and the 
resulting expression is: 


x(t) = A cos y - hA® (3y sin y + y* cos y) + aE y sin y 
Buz we (64 


- *Y cos y + ay sin y + Yor COS y| ~ eee 


(303906) 
where, for convenience, y = wot. Another approximate solution is given 


by Cunningham [ref. 10], using the Ritz method. The result is: 


Kire= A cos \fug + 0.75hA2 ¢ (me . 17) 


Using the arbitrary values w= 1, h = 1, and A = 1, a comparison of 
these two solutions is shown in Table III and a graph in Figure 3. 
With the arbitrary values selected, the new method yields good results, 
when compared to the Ritz method, up to about one second. Examination 
of the series shows that better results could be obtained over a longer 
period of time if the quantity A were made less than unity. 

In the method originally proposed by Baycura [ref. 1], the trans- 
formed equation was the same as that which would result if the relation 


X(s) equals A/s were used again in Equation (3.3.2): 


X(s) = =a (3.3.18) 


For this expression, then the inverse transform is: 

mie) = cos \ie. t (3.3.19) 
assuming that A = h = Wy = 1. At t = 2, Ki) = 202951, andmat t = 1, 
x(t) = 0.156. Thus, if desired, the equation in s can be reduced to a 
linear expression, and this will yield a solution which is not in the 


form of an infinite series. 
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TABLE IIT 


Ritz Laplace 
Time Method Method Errow 
0.0 1.000 1.000 0.000 
0.2 0.965 0.960 0.005 
0.4 0.863 0.844 0.019 
0.6 0.701 0.658 0.043 
0.8 0.490 0.416 0.074 
pu’ 0.245 0.134 02a 
lec -0.017 -0.173 0.156 
1.4 -0.278 -0.485 0.207 
6 -0.519 -0.786 0.267 
[wsse: -0.724 -1.060 0.336 
250 -0.880 -1.290 0.410 


3.4 An Equation with a Derivative Raised to a Power 


Consider now an equation of the form [ref. 11]: 
x(t) + [x(t)]2 + x(t) + k = 0, x(0) = 1, x(0) = 0 (oral) 


The transformed equation is found as before, thus: 
s*X2(s) - (s?-1)X(s) +£ = 0 (3.4.2) 


This is then solved for X(s) as before, and the result is: 


X(s) = a + \ - ates J (3.4.3) 


26 


The term in brackets is then expanded as a binomial series, but this 
time the negative sign in front of the radical is chosen to avoid a 


delta function in the solution, so that: 


sek Es? chase 
he aD kis G5) Ela (3.4.4) 


The inverse transform of the first term is k sinh t, but for the suc- 
ceeding terms the same difficulty is encountered as in Section 3.3. 


However, the table of transforms [ref. 12] lists: 


 — ven ty F(t) 
(C0 ee (3.4.5) 


where 1, a(t) is the modified Bessel function [ref. 13] given by: 
Z 


N(t)  & KEP (N+K+T) ial 


For the second term of the series, n = 3, N = 3/2, and: 


a 
I3(t) = ) (e/a) 
Z k=o St 
3 ae AL 
2 2 2 
- Myayr , Lefay* (e/a) Aon (3.4.7) 


and to evaluate the Gamma function [ref. 14]: 


a) ee (3.4.8) 

Thus: e a a 
2 Z 

ra(t) = 4t/2)* , B(t/2)" , 16(t/2) * , (3.4.9) 

z 3VT 15 Vr 105y1 
This result is substituted in Equation (3.4.5) to get: 

-] S wel Ot oat eas ie * 

Da fear = 35 | oe Toe 7 (3.4.10) 


Examination of this expression shows that the system is unstable, since 
the terms are always positive and the series has no limit, increasing 


without bound as t becomes large. 
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CHAPTER 4 


4.0 Conclusions 

It has been established that a method for obtaining an approximate 
solution to a non-linear differential equation has been derived, and 
that the method uses techniques and characteristics inherent to the 
Laplace Transform. Further, it has been shown that sufficient accuracy 
exists in the approximation to make it useful in many engineering 
applications. 

As in any approximation which uses a truncated infinite series to 
obtain a numerical answer to a problem, extreme care must be exercised 
to ensure that the approximation is accurate within a reasonable degree. 
This means that the parameters must be examined closely to avoid 
exceeding the capability of the approximation. It also means that the 
approximation may indeed only be accurate over a small part of the 
entire solution. 

The results included here are the outcome of what really amounts to 
preliminary investigation only. One question that remains is that, 
Since the solutions obtained in Sections 3.1 and 3.2 are so similar in 
form to the exact solutions, and since the solution seems exact in the 
Laplace expression, why isn't it, in fact, exact? After all, one major 
attraction of the Laplace Transform in solving linear equations is that 
the solution may be obtained by manipulating the transformed equation 
algebraically, and that is how the solutions were arrived at here. Thus, 
it seems that further research must be done in the area of this non- 


linear transform. 


28 


BIBLIOGRAPHY 


Baycura, 0. M., Approximate Laplace Transforms For Non-linear 
Differential Equations, Asilomar Conference on Circuits and 
Systems, Proceedings, Pacific Grove, California, 1968. 

Baycura, 0. M., Laplace Transforms for Quadratic Non-linearities, 
Submitted to, Journal of Mathematical Analysis and Applications, 
University of California, Berkeley. 

Roberts, G. E. and H. Kaufman, Table of Laplace Transforms, p. 223, 
equation (18), W. B. Saunders Co., Philadel chia. 1966. 

Cunningham, W. J., Introduction to Non-linear Analysis, p. 124, 
McGraw-Hill, New York, 1958. 

fond, 6. 163. 


Heerner, S. F., Fluid-Dynamic Drage P.cliqSomh@ennehyeSmediee, 
Midland Park, N. J., 1958. 


Roberts, G. E. and H. Kaufman, op. cit., p. 213, equation (67). 
Ibid., p. 170, equation (12). 


th 


Selby, S. M., Standard Mathematical Tables, 15° ed., p. 464, 


Chemical Rubber Co., Cleveland, 
Cunningham, W. J., op. cit., p. 161. 
Thaler, G. J., and M. P. Pastell, Analysis and Design of Non-linear 

Feedback Control Systems, p. 13/7, problem 3-3e, McGraw-Hill, 1962. 
Roberts, G. J., and H. Kaufman, op. cit., p. 213, equation (68). 


Selby, Stem. , op.-eit., p. 465, #iom 
IbiGiag~pe393 4478. 


29 


INITIAL DISTRIBUTION LIST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library 
Naval Postgraduate School 
Monterey, California 93940 


Naval Electronics and Communications Systems Command 
Washington, D. C. 


Professor 0. M. Baycura 

Department of Electrical Engineering 
Naval Postgraduate School 

Monterey, California 93940 


Lieutenant Charles R. Brady, USN 
335 First Street 
Palisades Park, New Jersey 07650 


Lieutenant Claus E. Zimmerman 
Box 1870 

Naval Postgraduate School 
Monterey, California 93940 





30 


No. Copies 


20 


UNCLASSIFIED 


Security Classification 


DOCUMENT CONTROL DATA-R&D 


(Security classification of title, body of ebetract and indexing ennotation muat be entered when the overall report ia classified 
1. ORIGINATING ACTIVITY (Corporate author) 2a. REPORT SECURITY CLASSIFICATION 


UNCLASSIFIED 


25. GROUP 


Naval Postgraduate School 
Monterey, Californta 93940 


3 REPORT TITLE 


Approximate Solutions to Non-Linear Differential Equations Using 
Laplace Transform Techniques 


' 4. OESCRIPTIVE NOTES (Type of report and, inclueive dates) 


Master's Thesis 


8. AUTHOR(S) (Firet name, middie initial, leet name) 


Charles R. Brady, Lieutenant, USN 


7. REPORT DATE 7a. TOTAL NO. OF PAGES 7b. NO. OF REFS 
“hprtt 1969 


Be. CON TRACT OR GRANT NO. Re. OCRIGINATOR’S REPORT NUMBER(S) 










—((. PROJECT NO. 





9b. OTHER REPORT NO(S) (Any other numbers that may be assigned 


c. 
thie report) 









d. 
/10. OISTRIGUTION STATEMENT 
This document has been approved for public release and sale; its distribution 
is unlimtted. 










12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 
Monterey, California 93940 





11. SUPPLEMENTARY NOTES 








| 13. A@STRACT 


At the present time the primary method of obtaining solutions to non-linear 
differential equations is by means of the digital computer and numerical techniques. 
A method 1s here proposed to find an approximate mathematical expression through the 
use of Laplace Transform techniques. Thus, the Laplace Transform concept 1s 


extended to the solution of non-linear differential equations. 


DD °"..1473 (Pace 1) 
S/N 0101-807-6811 3] ch CTassilicetion 


4- 31408 


UNCLASSIFIED 


Security Classification 





: Snore] of MESS | noc elle 
Non-linear Laplace Transforms 
Approximate solution 
Non-linear differential equations 
: 
: 
FORM 
DD "..1473 (8ack) UNCLASSIFIED | 
et 6 IN Security Classification a 




















































































































i 













































































thesB7972 
pproximate solutions to non-linear diff 


DUDLEY KNOX LIBRARY 


